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Abstract: 
spacetime evolution of quark-gluon plasma (QGP). The behavior of the speed of sound with respect to temperature and 


The speed of sound in quark matter 1s an important physical quantity for studying the properties and the 


density can reveal to some extent the equation of state and the phase structure of QGP. Building upon the previous studies 
on the speed of sound in symmetric quark matter, the formulae for calculating the speed of sound in asymmetric quark 
matter in the temperature-density space are further derived. The PNJL model is then used to calculate the dependence 
of the speed of sound on isospin asymmetry. Furthermore, the relationship between the magnitude of the speed of sound 
and the QCD phase structure 1s discussed, and the regions where the acoustic equation fails are indicated under different 


physical conditions. 
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Introduction 


Exploring the equation of state (EOS) and the phase 
structure of QGP is a hot topic in nuclear physics!-?), The 
speed of sound in QGP is a crucial mechanical characteristic 
to describe the variation of EOS. Its dependence on tempera- 
ture and baryon density carries important information in de- 
scribing the evolution of a fireball and the final observables. 
Recently, there have been some attempts to extract the speed 
of sound of QGP from the experimental data. For example, 
the studies indicate that the speed of sound is a function of 
charged particle multiplicity (dN,,,/dn)"°'!. An interest- 
ing method to estimate the speed of sound is proposed by 
building its connection with the net baryon number cumu- 
lants and the QCD phase structure!®!. 

Besides, there has been growing interest in the speed of 


[921] The variation of the speed 


sound in neutron star matter 
of sound with baryon density affects the mass-radius relation 


and the tidal deformability of neutron stars, and provides a 
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useful probe for studying the EOS of neutron star matter. To 
obtain a neutron star with a mass greater than twice that of the 
sun, some studies show that the EOS of neutron star matter 
must be sufficiently stiff in some density range, correspond- 
ing to a speed of sound squared much greater than 1/3 ??-281, 
The speed of sound in neutron star matter is also connected 
to the gravitational wave frequencies induced by the g-mode 


oscillation of a neutron starU?l, 


Furthermore, the magni- 
tude of speed of sound affects the dynamics of primordial 
density perturbations in the era of cosmic QCD phase tran- 
sition, which is relevant for the induced gravitational wave 
signals P’, 

Some calculations about the speed of sound in 
QCD matter have been performed in lattice QCD?!*!, 
(P)NJL model P640], TERT 


hadron resonance gas (HRG) model! 


quark-meson coupling model! 
42-43] field correlator 
method (FCM) 4*5] and quasiparticle model °!, In particu- 
lar, the relationship between the speed of sound and the QCD 
phase transition has been intensively studied ^5!, The nu- 
merical results indicate that the dependence of the speed of 
sound on the system parameters (temperature, density and 
chemical potential) is closely related to the phase structure 
of QCD. 

By far, the studies in the literature mainly focus on the 


speed of sound in symmetric quark matter, and there is a lack 


qi 
dq 
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of systematic research on the speed of sound in the isospin 
asymmetric quark matter. In this work we will explore the 
behavior of the speed of sound in asymmetric quark matter. 
On the one hand, we will derive the formulae for calculat- 
ing the speed of sound in asymmetric quark matter with two 
conserved charges. On the other hand, we will demonstrate 
numerically the dependence of the speed of sound on isospin 
asymmetric parameter in the PNJL model. The relevant re- 
search has certain significance in exploring the equation of 
state of QCD and the evolution of QGP. 


1 Formula of speed of sound in asymmetric 


quark matter and the PNJL quark model 


The general definition of speed of sound is ve 三 ( 也 P 
A specifying constant quantity X 1s required to describe the 
propagation of the compression wave through a medium. 
For a fireball created in relativistic heavy-ion collisions, it 
evolves with a constant entropy density to baryon density ra- 
tio s/pp under the ideal fluid approximation. Therefore, the 


speed of sound can be calculated along the isentropic curve 


0 
a t (1) 
sip g ðE A 


In practice the above definition can only be used to cal- 
culate the speed of sound along isentropic trajectories. To 
calculate the speed of sound throughout the full phase dia- 
gram, it is necessary to derive the corresponding formula as 
a function of temperature and baryon density. With the fun- 
damental thermodynamic relation, Eq. (1) can be rewritten 


as 
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where o = s/pp is defined for convenience. 


_ 0p, 6) 
E (T. at 


The speed of sound will be calculated under the physical 
condition p, — 0 for strange quark, which is the requirement 
of strangeness conservation in strong interaction. The asym- 
metric parameter a is defined as a = p,/p, to describe the 
ratio of u, d quark number density. Then, we can derive the 


following formulae in the grand canonical ensemble 
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For a given asymmetric parameter a, Eq. (2) can be used to 
calculate the speed of sound in the T — pp plane with the 
above formulae. For the case of a = 1, these formulae return 
to those for symmetric quark matter. 

In the numerical calculation, we take the 2+1 flavor 
PNJL quark model. the Lagrangian density is written as 


8 


L = Giy"D,+7oh-im)q+G Y, [Gag iy, A0] 
k=0 
-K |det ,(G(1 + y5)q) + det (GU. — 75)q)| 


where q denotes the quark fields with three flavors, u, d, and 
5; fiy = diag(m,, mg, m,) in flavor space; G and K are the 
four-point and six-point interacting constants, respectively. 
The f = diag(u,. Wa, Hs) are the quark chemical poten- 
tials. The covariant derivative in the Lagrangian is defined 
as D, = 0, —iA,. The gluon background field A, = 89 Ag is 
supposed to be homogeneous and static, with Ay = g.A5 Z ; 
where £ is the SU (3) color generators. The effective poten- 


tial U(GIA ], ®[ A], T) is expressed with the traced Polyakov 


loop ® = (Tr, L)/ Nç and its conjugate ® = (Tr, LY Nc. 
The Polyakov loop L is a matrix in color space 
B 
By = Pexp i | dc Ak, 可 (12) 


0 
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where p = 1/T is the inverse of temperature and A, = i Ag. 


The Polyakov-loop effective potential is 


- - Dao + kryn[1 - 660 (13) 


+4(B? + D°) - 360, 


where 

a(T)- aga, (2) sa (2) and b(T)=b, ey (14) 
The parameters ay = 3.51, a, = —247, a; = 15.2, b = 
—1.75 and Ty = 270 MeV are fitted according to the earlier 
lattice simulation of QCD thermodynamics in the pure gauge 
sector '?!, When the quark fields are included T, = 210 MeV 
will be implemented to describe well the thermodynamics 
and phase transition of QCD. The constituent quark mass in 


the mean field approximation can be derived as 
Mi =m, —4G¢q; + 2K od, (i zx jk) (15) 


where d$; stands for quark condensate of the flavor i. 


The thermodynamic potential of bulk quark matter is de- 


rived as 
Q = U($,0,T)-2G (6 + oy’ + 6) - 4K$,9, 40, 
-2| SPE 十 五 十 五 ) 
(2? u d s 
A 
dp 
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2) ay * 9) (16) 
where Q) = ln( + 36e (FT + 3e EMT + 


p Mc HS Q, - In(1 十 3de (ErtuT 十 3@e XErtu)T 4 
e XErtWT). and E, = 


Mi is the chemical potential for the quark flavor i. The pres- 


4/ P? + M? is the dispersion relation. 


sure and energy density can be derived using the thermody- 


namic relations in the grand canonical ensemble as 


p=-Q, e=-p+Ts+ >, Hj pj (17) 


i-ud,s 

For the given temperature T', baryon density pp and 

asymmetric parameter a, the values of $,, pa, Ps P, ®, p, 

H4 and p, can be determined by solving the equations de- 
rived by minimizing the thermodynamic potential 

NQ 00 00 00 00 _ 

0$, dpa d$, d 8B 


In the numerical calculation, a cut-off A is implemented in 


0. (18) 


the three-momentum space for divergent integrations. The 
model parameters are taken with: A = 602.3 MeV, GA? = 
1.835, KA? = 12.36, M,a = 5.5 and m, = 140.7 MeV, 


determined by fitting f, = 92.4 MeV, M, = 135.0 MeV, 
my = 497.7 MeV and m, = 957.8 MeV P? , 


2 Numerical results and discussions 


We first present in Fig. 1 the isentropic curves for o = 
50, 10, 5, 1 in the T — ug plane. Each curve in this figure 
corresponds to an evolutionary trajectory of QGP generated 
with a specific collision energy. For the two cases of o = 50 
and 10, the fireball only passes through the stable phase dur- 
ing the evolution. However for the cases ofo = 5 and 1, both 
the metastable and unstable phases are involved. The behav- 
iors of o in the T — Ap plane is mainly determined by the 
phase structure of the system. In theory, for a fermion system 
with the temperature and density dependent dynamic mass, 
there usually exists a first-order phase transition at low tem- 
perature. Numerical calculations in different systems show 
that the qualitative behaviors of s/p g are similar to the curves 
shown in Fig. 1. For the liquid-gas phase transition in nu- 
clear matter, we also derived the similar curves of s/pg in 
Fig. 1(c) in Ref. 9 . Generally, the appearance of the inflec- 
tion point does not directly corresponds to the position ofthe 
phase transition, but the crossover phase transition of QCD 
at high temperature can induce some perturbation on s/p, as 
shown in the curve of o = 50. 
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Fig. 1 (Color online) Isentropic curves for o. = 50,10,5,1 for 


both the symmetric and asymmetric quark matter in the 
T — ug plane. 


We plot in Fig. 2 the EOSs of quark matter along differ- 
ent evolutionary paths given in Fig. 1. The numerical results 
show that for the cases of o = 1 and 5 in symmetric matter 


and o = 1 in asymmetric matter, the pressure of the system 
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has a maximum and a minimum as the energy density in- 
creases, respectively. The region between the maximum and 
minimum approximately correspond to (lies in exactly) the 
unstable phase ofthe system. After passing through the mini- 
mum, the pressure grows rapidly with the increase of energy 
density. For the cases of o = 5 in asymmetric quark mat- 
ter and o = 10 and 50 in both the symmetric and asymmet- 
ric quark matter, the pressure increases monotonically with 
the rising energy density. Fig. 2 also indicates that the equa- 
tions of state of asymmetric matter are stiffer than those of 
symmetric matter, although it is not obvious for the case of 
o — 50. These behaviors of the equations of state under dif- 
ferent physical conditions determine the speed of sound in 


quark matter. 
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Fig.2 (Coloronline) EOSs along the different evolutionary paths 
shown in Fig. 1. 


The evolving behaviors of the speed of sound corre- 
sponding to the EOSs under different conditions in Fig. 2 
are plotted in Fig. 3. This figure demonstrates that the in- 
fluence of isospin effect on the speed of sound of quark mat- 
ter. In the chiral restored phase (the high-temperature and/or 
high-density region), the isospin splitting effect is relatively 
smaller, because the differences of both the quark masses and 
chemical potentials between u and d quarks are very small, 
leading to that the asymmetric parameter a has a weak effect 
on the EOSs. The isospin effect is also weak at lower den- 
sity (lower energy density) in the chiral broken phase with 
the similar reasons for the chiral restored phase, in particular 
when the mass ofu quark is the same as that of d quark. How- 
ever, the isospin effect has a relatively larger influence on the 
EOSs at intermediate density and in the phase transition re- 


gion. Therefore, from the low-temperature and low-density 
region to the high-temperature and/or high-density region, 
the influence of the isospin effect on EOSs goes from weak 
to strong and then to weak again, which is consistent with the 
results shown in Fig. 2. With the combination of the above 
qualitative analysis and the EOSs in Fig. 2, as well as the 
Isetropic trajectories in Figs. 1, 4 and 5, it can be understood 
that when o is relatively smaller (e.g., o —1, 5, 10), the speed 
of sound of asymmetric quark is larger than that of symmet- 
ric quark matter in the phase transition region, and a reversal 
happens when approaching the chiral restored phase. while 
for o. = 50, the sound speed curves of the symmetric and 
asymmetric quark mater almost coincide because the isen- 
tropic curves of them mainly pass through the low-density 
region and the chiral restored phase at high temperature. 
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Fig. 3 (Color online) Square of speed of sound for the EOSs under 
different conditions along the evolutionary paths shown in 
Fig. 1. 


To demonstrate the dependence of the speed of sound on 
temperature and density, we present the contour plots of the 
speed of sound in both the symmetric quark matter and asym- 
metric quark matter in Figs. 4 and 5, respectively. From the 
two figures, it can be seen that in the low-temperature region 
of the first-order phase transition, the restoration of chiral 
symmetry causes the square of the speed of sound to quickly 
approach 1/3, as indicated by the red curve in Fig. 3. In the 
smooth transition region at high temperatures, the growth of 
the speed of sound with increasing temperature or density 
is relatively slower. However, at extremely high tempera- 
tures and densities, the square of the speed of sound still ap- 
proaches 1/3. 


The red solid lines in Figs. 4 and 5 correspond to the 
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boundaries where the speed of sound becomes zero with the 
condition of (2) = 0. As pointed out in our previous 
study 7, the acoustic wave equation becomes a decaying 
function within the region enclosed by the red solid lines, 
and the disturbances cannot propagate outwards within this 
region. Comparing Fig. 4 and 5, we observe that the region 
with the negative speed of sound squared in asymmetric mat- 
ter is smaller than that in symmetric matter, which can also 


be seen from Fig. 3. 
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Fig. 4 (Color online) Contour map of c? in the T — pp plane for 


symmetric quark matter with a = 1. The red curve is the 


boundary of vanishing sound speed. 
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Fig. 5 (Color online) Contour map of c? in the T — p, plane for 
asymmetric quark matter with a = 1/2. The red curve is 
the boundary of vanishing sound speed. 


Furthermore, combined with the contour plots of the 
speed of sound under the condition p, = py = ju, 7l, it can 
be observed that the condition p, = 0 taken in this study has 
a significant influence on the speed of sound in quark mat- 
ter, especially in the low temperature and high density region 
where the chiral symmetry of u and d quarks has restored. 
The calculation" demonstrates that under the condition of 


H, = Hg = Mg, the speed of sound at low temperatures rapidly 
increases to nearly TE after the restoration of chiral sym- 
metry of u and d quarks. However, as the density increases, 
the presence of strange quarks causes the equation of state 
to become softer, resulting in a subsequent decrease of the 
sound speed, which tends to be V13 again after the restora- 
tion of chiral symmetry of s quark at higher densities. 

To compare the value of the speed of sound in sym- 
metric and asymmetric quark matter, we plotted the contour 
maps ofr, = E p! a in Fig. 6. In the gray area, the speed 
of sound in asymmetric quark matter is slightly larger than 
that in symmetric quark matter, andr, = 1 on the boundary of 
the gray area. In the red region, r, < 1, it means that the value 
of the speed of sound in asymmetric quark matter is smaller 
than that in symmetric quark matter. It should be noted that 
the difference in speed of sound between the two cases in 
the gray area is very small (r, < 1.005 in most of the gray 
region). The negative ratio appears in the region between 
the two boundaries (black solid curve and dashed curve) of 
the vanishing sound speed. In this region, the square of the 
speed of sound in asymmetric quark matter is positive, while 
in symmetric quark matter, the mechanical instability vio- 
lates the sound wave equation with the c? being negative. 
The region where the ratio takes a larger absolute value lo- 
cates near the boundary of (black solid curve) of the vanish- 
ing sound speed in symmetric quark matter. 
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Fig. 6 (Color online) Contour map of r, = pes wl Ba in the T — 


Pp plane. The solid (dashed) black curve is the boundary of 
vanishing sound speed in symmetric (asymmetric) matter. 
In the gray area, the squared speed of sound in asymmetric 
matter is slightly larger than that in the symmetric matter 
with r, > 1;r, = 1 on the boundary ofthe gray area; r, < 1 
in the red area. 
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3 Summary 


We derived the formula for calculating the speed of 
sound in asymmetric quark matter in the temperature and 
density space, and calculated the speed of sound in the PNJL 
quark model. The results indicate that the boundary of van- 
ishing sound speed in asymmetric quark matter is smaller 
than that in symmetric quark matter, meaning that the range 
where the acoustic wave equation fails in asymmetric quark 
matter is smaller than that in symmetric quark matter. The 
calculations also indicate that in most of the stable phase, the 
speed of sound in asymmetric quark matter is slightly larger 
than that in symmetric quark matter. Additionally, it should 
be pointed out that only the change of sound speed caused by 
isospin asymmetry of u, d quark chemical potential is con- 
sidered in this study. The difference in u, d quark dynamic 
masses is not included. After introducing the isospin inter- 
actions that cause the splitting of u, d quark masses, it is ex- 
pected that the isospin effects on the speed of sound will be 
more pronounced. The relevant research is currently under- 


way. 


4 Acknowledgments 


The work was supported by National Natural Science 
Foundation of China (Grant No. 11875213) and Natural Sci- 
ence Basic Research Plan in Shaanxi Province of China (Pro- 
gram No. 2024JC-YBMS-018). 


References 


1] FUKUSHIMA K. Phys Rev D, 2008, 77: 114028. DOI: 10.1103/Ph 
ysRevD.77.114028. 

2] RATTIC, THALER M A, WEISE W. Phys Rev D, 2006, 73: 014019. 
DOI: 10.1103/PhysRevD.73.014019. 

3] COSTA P, RUIVO M C, DE SOUSA C A, et al. Symmetry, 2010, 2 
(3): 1338. DOI: 10.3390/sym203 1338. 

4] FU WJ, ZHANG Z, LIU Y X. Phys Rev D, 2008, 77: 014006. DOI: 
10.1103/PhysRevD.77.014006. 

5] SASAKI T, TAKAHASHI J, SAKAI Y, et al. Phys Rev D, 2012, 85: 
056009. DOT: 10.1103/PhysRevD.85.056009. 

6] FERREIRA M, COSTA P, PROVIDÉNCIA C M C. Phys Rev D, 
2014, 89: 036006. DOI: 10.1103/PhysRevD.89.036006. 

7] SCHAEFER B J, WAGNER M, WAMBACH J. Phys Rev D, 2010, 
81: 074013. DOI: 10.1103/PhysRevD.81.074013. 

8] SKOKOV V, FRIMAN B, REDLICH K. Phys Rev C, 2011, 83: 
054904. DOI: 10.1103/PhysRevC.83.054904. 

9] QIN S X, CHANG L, CHEN H, et al. Phys Rev Lett, 2011, 106: 
172301. DOI: 10.1103/PhysRevLett.106.172301. 

GAO F, CHEN J, LIU Y X, et al. Phys Rev D, 2016, 93: 094019. 
DOI: 10.1103/PhysRevD.93.094019. 


[ 


11 


[12 


[13 


[14 


[15 


[16 


[17 


18 


[19 


[20 


21 


[22 


[23 


[24 


[28 


[29 


[30 


[31 


[32 


[33 


34 


35 


36 


FISCHER C S, LUECKER J, WELZBACHER C A. Phys Rev D, 
2014, 90: 034022. DOI: 10.1103/PhysRevD.90.034022. 

FU W J, PAWLOWSKI J M, RENNECKE F. Phys Rev D, 2020, 101: 
054032. DOI: 10.1103/PhysRevD.101.054032. 

WANG X, WEI M, LI Z, et al. Phys Rev D, 2019, 99: 016018. DOT: 
10.1103/PhysRevD.99.016018. 

LIU L M, XU J, PENG G X. Phys Rev D, 2021, 104: 076009. DOI: 
10.1103/PhysRevD.104.076009. 

WEN X J, HE R, LIU J B. Phys Rev D, 2021, 103: 094020. DOI: 
10.1103/PhysRevD.103.094020. 

GARDIM F G, GIACALONE G, LUZUM M, et al. Nature Physics, 
2020, 16(6): 615. DOI: 10.1038/s41567-020-0846-4. 

BISWAS D, DEKA K, JAISWAL A, et al. Phys Rev C, 2020, 102: 
014912. DOT: 10.1103/PhysRevC.102.014912. 

SORENSEN A, OLIINYCHENKO D, KOCH V, et al. Phys Rev Lett, 
2021, 127: 042303. DOI: 10.1103/PhysRevLett. 127.042303. 

REED B, HOROWITZ C J. Phys Rev C, 2020, 101: 045803. DOI: 
10.1103/PhysRevC.101.045803. 

KANAKIS-PEGIOS A, KOLIOGIANNIS P S, MOUSTAKIDIS C C. 
Phys Rev C, 2020, 102: 055801. DOI: 10.1103/PhysRevC.102.055 
801. 

HAN S, PRAKASH M. The Astrophysical Journal, 2020, 899(2): 
164. DOI: 10.3847/1538-4357/aba3c7. 

TEWS I, CARLSON J, GANDOLFI S, et al. The Astrophysical Jour- 
nal, 2018, 860(2): 149. DOI: 10.3847/1538-4357/aac267. 

GREIF S K, RAAIJMAKERS G, HEBELER K, et al. Monthly No- 
tices of the Royal Astronomical Society, 2019, 485(4): 5363. DOI: 
10.1093/mnras/stz654. 

FORBES M M, BOSE S, REDDY S, et al. Phys Rev D, 2019, 100: 
083010. DOI: 10.1103/PhysRevD.100.083010. 

DRISCHLER C, HAN S, LATTIMER J M, et al. Phys Rev C, 2021, 
103: 045808. DOI: 10.1103/PhysRevC.103.045808. 

ESSICK R, TEWS I, LANDRY P, et al. Phys Rev C, 2020, 102: 
055803. DOI: 10.1103/PhysRevC.102.055803. 

HAN S, MAMUN M A A, LALIT S, et al. Phys Rev D, 2019, 100: 
103022. DOI: 10.1103/PhysRevD.100.103022. 

KOJO T. AAPPS Bulletin, 2021, 31(1): 11. DOI: 10.1007/s43673-0 
21-00011-6. 

JAIKUMAR P, SEMPOSKI A, PRAKASH M, et al. Phys Rev D, 
2021, 103: 123009. DOI: 10.1103/PhysRevD.103.123009. 

ABEK T, TADA Y, UEDA I. Journal of Cosmology and Astroparticle 
Physics, 2021, 2021(06): 048. DOI: 10.1088/1475-7516/2021/06/04 
8. 

AOKI Y, ENDRGÓDI G, FODOR Z, et al. Nature, 2006, 443(7112): 
675. DOI: 10.1038/nature05120. 

BORSANYIS, FODOR Z, HOELBLING C, et al. Phys Lett B, 2014, 
730: 99. DOI: https://doi.org/10.1016/j.physletb.2014.01.007. 
BAZAVOV A, BHATTACHARYA T, DETAR C, et al. Phys Rev D, 
2014, 90: 094503. DOI: 10.1103/PhysRevD.90.094503. 
PHILIPSEN O. Progress in Particle and Nuclear Physics, 2013, 70: 
55. DOI: https://doi.org/10.1016/j.ppnp.2012.09.003. 

BORSANYI S, FODOR Z, GUENTHER J N, et al. Phys Rev Lett, 
2020, 125: 052001. DOI: 10.1103/PhysRevLett.125.052001. 
MOTTA M, STIELE R, ALBERICO W M, et al. The European Phys- 
ical Journal C, 2020, 80(8): 770. DOI: 10.1140/epjc/s10052-020-82 
18-x. 


原子 核 物 理 评 论 


第 0 卷 


GHOSH S K, MUKHERJEE T K, MUSTAFA M G, et al. Phys Rev 
D, 2006, 73: 114007. DOI: 10.1103/PhysRevD.73.114007. 
MARTY R, BRATKOVSKAYA E, CASSING W, et al. Phys Rev C, 
2013, 88: 045204. DOI: 10.1103/PhysRevC.88.045204. 

DEB P, KADAM G P, MISHRA H. Phys Rev D, 2016, 94: 094002. 
DOI: 10.1103/PhysRevD.94.094002. 

ZHAO Y P. Phys Rev D, 2020, 101: 096006. DOI: 10.1103/PhysRe 
vD.101.096006. 

ABHISHEK A, MISHRA H, GHOSH S. Phys Rev D, 2018, 97: 
014005. DOT: 10.1103/PhysRevD.97.014005. 

VENUGOPALAN R, PRAKASH M. Nuclear Physics A, 1992, 546 
(4): 718. DOI: https://doi.org/10.1016/0375-9474(92)90005-5. 
BLUHM M, ALBA P, ALBERICO W, et al. Nuclear Physics A, 2014, 
929: 157. DOI: https://doi.org/10.1016/j.nuclphysa.2014.06.013. 
KHAIDUKOV Z V, LUKASHOV M S, SIMONOV Y A. Phys Rev 
D, 2018, 98: 074031. DOI: 10.1103/PhysRevD.98.074031. 
KHAIDUKOV Z V, SIMONOV Y A. Phys Rev D, 2019, 100: 
076009. DOT: 10.1103/PhysRevD.100.076009. 

MYKHAYLOVA V, SASAKI C. Phys Rev D, 2021, 103: 014007. 
DOI: 10.1103/PhysRevD.103.014007. 

HE WB, SHAO GY, GAO X Y, et al. Phys Rev D, 2022, 105: 094024. 
DOI: 10.1103/PhysRevD.105.094024. 

HE W B, SHAO GY, XIE C L. Phys Rev C, 2023, 107: 014903. DOL 
10.1103/PhysRevC.107.014903. 

ROBNER S, RATTI C, WEISE W. Phys Rev D, 2007, 75: 034007. 
DOI: 10.1103/PhysRevD.75.034007. 

REHBERG P, KLEVANSKY S P, HÜFNER J. Phys Rev C, 1996, 53: 
410. DOI: 10.1103/PhysRevC.53.410. 


第 0 期 XIE Chonglong et al: Speed of sound in asymmetric quark matter... 


非 对 称 夸 元 物质 中 的 声速 


AER, HRR, SSH, 8p 


inl 


H 


(1. 西安 交通 大 学 ， 


12x 710049; 


2. 北京 大 学 ， 北 京 100871 ) 


in 


摘要 : 夸克 物质 中 的 声速 是 研究 压 克 和 胶 子 等 离子 体 的 性 质 和 时 空 演化 的 一 个 重要 物理 量 ， 声 速 随 温 度 和 密度 
的 变化 行为 在 一 定 程度 上 能 够 揭示 QGP 的 物 态 和 相 结构 。 在 之 前 研究 对 称 夸克 物质 中 声速 的 基础 上 ， 本 文 进 


一 步 推 导 了 在 温度 和 密度 空间 计算 非 对 称 夸克 物质 声速 的 公式 ,六 
对 同位 旋 不 对 称 参数 的 依赖 ,讨论 了 声速 的 大 小 与 QCD 相 结 构 的 关系 ,指出 了 不 同 条 件 下 声波 方程 失效 的 范 


S 


Hl 。 
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F 利 用 PNJL 模型 计算 了 夸克 物质 中 的 声速 
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